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Assistant Commissioner for Patents, 
Washington D.C. 20231 



Sir: 



Applicants herewith amend the above-referenced PCT application, and 
request entry of the Amendment prior to examination on the United States 
Examination Phase. 



On page 21: 

replace line 1 with -WHAT IS CLAIMED IS:-; 

Please replace original claims 1-10 with the following rewritten claims 1-10, 
referring to the mark-ups in Appendix A. 

1 . (Amended) A method for designing a technical system, comprising the 
steps of: 

a) providing a substitute model that describes measurement data of a 
predetermined system; 

b) determining a numerical value for a quality of said substitute model by 
comparing said measurement data of said predetermined system with data 
determined by said substitute model; 



IN THE CLAIMS: 
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c) adapting said substitute model from said numerical value for said quality to 
be as higli of a quality as possible; 

d) applying said substitute model adapted with regard to its quality in a design 
of said technical system. 

5 

2. (Amended) The method as claimed in claim 1, wherein said substitute 
model is a regression model. 

3. (Amended) The method as claimed in claim 1 , wherein said step of 

10 determining a numerical value for a quality further utilizes a mean square deviation 
of said measurement data from said data determined by said substitute model. 



4. (Amended) The method as claimed in claim 1 , further comprising the step 

of: 

15 sorting said measurement data according to their quality, with respect to the 

deviation of the latter from said data determined by said substitute model; and 
picking out a predetermined number of n% of worst measurement data. 

5. (Amended) The method as claimed in claim 1 , further comprising the step 

20 of: 

sorting said measurement data according to their quality, with respect to the 
deviation of the latter from said data determined by said substitute model; and 

picking out a predetermined number of n% of worst measurement data unless 
this data lie in a continuous range. 

25 

6. (Amended) The method as claimed in claim 1 , further comprising the step 

of: 

reducing an amount of measurement data in the course of a preprocessing 
operation. 

30 

7. (Amended) The method as claimed in claim 6, further comprising the step 

of: 
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classifying, in which said preprocessing operation, of said measurement data. 

8. (Amended) The method as claimed in claim 1 , further comprising the step 

of: 

controlling a technical plant utilizing said data obtained by designing. 

9. (Amended) The method as claimed in claim 8, further comprising the step 

of: 

online adapting control for said technical plant. 

10. (Amended) An arrangement for designing a technical system, 
comprising: 

a processor unit which is set up in such a way that 

a) measurement data of a predetermined system are described based on a 
substitute model and stored in said processor unit; 

b) a numerical value for a quality of said substitute model is determined by 
said processor unit by comparing said measurement data of the predetermined 
system with data determined by said substitute model; and 

c) said substitute model is adapted, utilizing said processor unit, from said 
numerical value for said quality to be as of high a quality as possible, wherein said 
substitute model adapted with regard to its quality is used for designing said 
technical system. 

REMARKS 

The present Amendment revises the specification and claims to conform to 
United States patent practice, before examination of the present PCT application in 
the United States National Examination Phase. Pursuant to 37 CFR 1 .125 (b), 
applicants have concurrently submitted a substitute specification, excluding the 
claims, and provided a marked-up copy. All of the changes are editorial and 
applicant believes no new matter is added thereby. The amendment, addition, 
and/or cancellation of claims is not intended to be a surrender of any of the subject 
matter of those claims. 
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Early examination on tlie merits is respectfully requested. 
Submitted by, 



Mark Bergner 
Schiff Hardin & Waite 
Patent Department 
6600 Sears Tower 
233 South Wacker Drive 
Chicago, Illinois 60606-6473 
(312)258-5779 
Attorneys for Applicant 
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Appendix A 
Mark Ups for Claim Amendments 

This redlined draft, generated by CompareRite (TM) - The Instant Redliner, shows the differences between - 
5 original document : Q :\DOCUMENTS\YEAR 200 1 \P0 1 0 1 58-SCHAEFFLER-DESIGNING A TECHNICAL 
SYSTEM\ORIGINAL CLAIMS.DOC 

and revised document: Q:\DOCUMENTS\YEAR 2001\P010158-SCHAEFFLER-DESIGNING A TECHNICAL 
SYSTEMVAMENDED CLAIMS.DOC 

10 CompareRite fomid 71 change(s) in the text 

Deletions appear as Overstrike text siurounded by [] 
Additions appear as Bold-Underline text 

15 1 . (Amended) A method for designing a technical system, comprising the 

steps of: 

r a) providin g a substitute model that describes [a) in which] measurement 

data of a predetermined system [ar e d e scribed on tho basis of a cubstituto modo l ;]i 
[b) in wh i ch] b) determining a numerical value for fthe} a quality of {the} said 
-20 substitute model [is dotorminod] by comparing {the} said measurement data of {the} 
said predetermined system with data determined by {the} said substitute model; 

c) [in which th o ] adapting said substitute model [is adapted] from [tho] said 
numerical value for {the} said quality to be fe^ as high of a quality as possible; 
2 d) [in which tho] applying said substitute model adapted with regard to its 

-25 quality [is used for doo i gning tho] in a design of said technical system. 

2. (Amended) The method as claimed in claim 1 , [in which tho] wherein 
said substitute model is a regression model. 

3. (Amended) The method as claimed in claim 1 [or 2, in which tho]^ . 
wlierein said step of determining a numerical value for a quality [is d o tormin o d 

30 on tho basis of] further utilizes a mean square deviation of {the} said measurement 
data from {the} said data determined by {the} said substitute model. 

4. (Amended) The method as claimed in [one of tho procoding c l aims, in 
which thol claim 1. further comprising the step of: 

sorting said measurement data [aro sorted] according to their quality, with 
35 respect to the deviation of the latter from {the} said data determined by {the} said 
substitute mode^i and 
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picking out a predetermined number of n% of {the} worst measurement data 
[arc picked out] . 

5. (Amended) The method as claimed in claim ['1, in wliich the n% of the ] r, 
further comprising the step of: 

sorting said measurement data according to their quality, with respect 
to the de viation of the latter from said data determined by said substitute 
model; and 

picking out a predete rmined number of n% of worst measurement data 
unless this data Faro not picknd nut 

i f they] lie in a continuous range. 

6. (Amended) Tiie method as claimed in [ono of the procoding claims, in 
which tho1 claim 1. further comprising the step of: 

reducing an amount of measurement data [is reduc e d] in the course of a 
preprocessing operation. 

7. (Amended) The method as claimed in claim 6, [in which tho ] further 
comprising the step of: 

classifyin g, in which said preprocessing operation [comprisos a 
classification of tho] , of said measurement data. 

8. (Amended) The method as claimed in [ono of the procoding claims, in 
which tho] claim 1. further comprising the step of: 

controlling a technical plant utilizing said data obtained by [moans of] 
designing [aro usod for controlling a tochnical plant.] . 

f& ^9. (Amended) The method as claimed in claim 8, ffo^ further 
comprising the step of: 

O'TiWne [adaptation of tho] adapting control for {the} said technical plant. 

1 0. (Amended) An arrangement for designing a technical system, [with] 
comprising: 

a processor unit which is set up in such a way that 

a) measurement data of a predetermined system are described based on 
[the basis of] a substitute model and stored in said processor unit : 
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b) a numerical value for {the] a quality of fthe} said substitute model is 
determined by said processor unit bv comparing fthe} said measurement data of 
the predetermined system with data determined by {the} said substitute model; and 

c) {the} said substitute model is adapted ffrom tho1 , utilizing said processor 
unit from said numerical value for {the} said quality to be as of {as} high a quality 
as possible , wherein said ff 

the} substitute model adapted with regard to its quality is used for designing 
{the} said technical system. 
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SPECIFICATION 

TITLE [DESCRIPTION] 

METHOD AND ARRANGEMENT FOR DESIGNING A TECHNICAL SYSTEM 

BACKGROUND OF THE INVENTION 

Field of the Invention 

1^°°°^^ The invention relates to a method and arrangement 

for designing a technical system. 

Description of the Related Art 

^0002] The system behavior of a technical system, for 

example_j^ a process engineering plant or system in heavy 
industry, depends on numerous parameters. In the course of 
designing such a system, [that ia in particular] particularly 
in the case of a new design or when adapting or adjusting an 
already existing system, it is necessary to comply with 
preconditions [, for example with regard to] such as the cost- 
effectiveness or environmental impact of the system. Each 
precondition is formulated as a target function, the 

[optimization of which with regard] general aim being to 
optimize this wit h respect to the other target functions 
the general aim. 

SUMMARY OF THE INVENTION 
I^QOQ^J The object of the invention is to make it possible 
for a technical system to be designed on the basis of 
measurement data of a predetermined system. Specifically with 
regard to the optimization of the existing system or with 
regard to an optimized new design of a system, such use of 
known measurement data is of great significance. 
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[0004] This object is achieved [according to the foaturca 
of the independent patent claimo . — Dcvclopmcnto of the 
invention aloo emerge from the dependent claims. 

method and apparatus described below. 

[0005] To achieve the object, the present invention 

provides a method for designing a technical system in which 
measurement data of a predetermined system are described on 
the basis of a substitute model [ia apecified] . A numerical 
value for the quality of the substitute model is determined by 
comparing the measurement data of the predetermined system 
with data determined by the substitute model. [On the basis 
ef-]- Based on the numerical value for the quality, the 
substitute model is adapted to be of as high a quality as 
possible . 

[0006] The substitute model adapted with regard to its 
quality is used for designing the technical system. 

[0007] The measurement data obtained from many different 
realized systems are used for describing the substitute model. 
[With the] The substitute model [it io attempted] attempts to 
replicate the predetermined system as well as possible. The 
numerical value for the quality of the replication is 
determined by comparing the actual measurement data with the 
data obtained on the basis of the substitute model. A great 
difference between the measurement data and the data of the 

substitute model corresponds to poor quality, that is_^ a poor 

mapping of the predetermined system into the substitute model. 
The numerical value for the quality is used to adapt the 
substitute model to make the quality itself become as high as 
possible and consequently to make the substitute model 
describe the predetermined system as well as possible. The 
high-quality substitute model obtained in this way is used for 
designing the technical system. 
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[QQ08] Designing is understood in a general sense as 
meaning both the new design of a technical system and the 
adaptation of an already existing technical system. 
[0009] One development comprises that the substitute model 

is a regression model. The regression model is based on the 
description 

[0010] Yi = fp(xi) + ei 

[0011] where 

[0012] (yi, Xi) denotes predetermined pairs of 

values , 

[0013] fp denotes a function which is dependent on a 

parameter p, and 

[0014] ei denotes an error. 
[0015] It is then intended to minimize the error (as a 

function of P) : 

Se? = <p(p). 

i = l 

[0016] 

[0017] If the following example 

[0018] Y = Po + pix + Psx^ + e 

[0019] is taken as a basis, the functional relationship is 

of a quadratic order, while the regression model (function, 
dependent on P) is linear. 

[0020] In another development, the quality can be determined 

on the basis of a mean square deviation of the measurement data 
from the data determined by the substitute model. The 
adaptation of the substitute model takes place by minimizing the 
mean square deviation. 
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[00213 One refinement comprises that the measurement data are 
sorted according to their quality, with respect to the deviation 
of the latter from the data determined by the substitute model , 
and a predetermined number of n% of the worst measurement data 
are picked out. Consequently, a quality is determined for each 
item of measurement data, the volume of measurement data, 
preferably in the form of a list, being sorted according to 
their quality and the n% of the worst measurement data, or the n 
worst measurement data, being picked out. In particular, it must 
be checked whether the n% of worst measurement data or n worst 
measurement data lie in a continuous range. If this is the 
case, these measurement data are not picked out, since they very 
probably do not represent measurement errors but a continuous 
range which is not mapped accurately enough by the substitute 
model . 

[0022] Another development [comrisca] comprises that the 
[mcaauro mcnt] measurement data are subjected to a 
preprocessing operation. Since, in an actual predetermined 
system, a large amount of measurement data occur per unit of 
time, it is advisable to subject these measurement data to a 
preprocessing operation and consequently ensure that largely 
significant measurement data are taken for forming the 
substitute model. The preprocessing operation preferably 
finds its expression in a reduction in the number of 
measurement data . 

[0023] In the process, the measurement data are divided 

into classes according to predetermined criteria. The 
measured values of a class are assessed and those measured 
values for which the assessment lies below a predetermined 
first threshold value are picked out. The picking out of the 
measured values results in a reduction with regard to the 
number of measured values. This results in a significantly 
reduced number of measured values for a further processing 
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operation. The further processing operation can take place 
with less computing effort in comparison with the unreduced 
number of measured values . 

The classes themselves can also be assessed. In 
particular, a class for which the assessment lies below a 
predetermined second threshold value can be picked out. As a 
result, an additional reduction of the number of measured 
values is obtained. 

10025} Another development of the preprocessing operation 

comprises that a criterion for the classification is that, for 
each class, measured values are determined as a default for 
setting parameters of the technical system. The technical 
system is typically set on the basis of a predetermined number 
of setting parameters; after setting, a reaction of the system 
to the setting parameters takes place (usually with a time 
delay) (transient response, transient phenomenon of the 
system) . After setting, consequently a certain set of 
measured values which can be assigned to the transient 
phenomenon are recorded, with measured values continuing to 
occur after the transient phenomenon has come to an end 

(transition to steady-state operation) , and are assigned to 
the predetermined set of setting parameters. [By adjuoting] 
Adjusting the setting parameters-b-+ determines a new class 
determined] . All of the measured values which respectively 
occur after adjustment of the setting parameters belong in a 
class of their own. 

1^0026] In addition, measured values of a class which can be 

assigned to the respective transient phenomenon can be picked 
out. Furthermore, erroneous measured values can be picked 
out. The setting of large technical systems is in many cases 
directed at long-term steady-state operation. It is advisable 
for the measured values which relate to the transient 
phenomenon (of short duration in relation to the steady-state 
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operation after the transient phenomenon has come to an end) 
to be picked out, since they have the effect of falsifying 
measured values for steady-state operation. In particular_j_ 
when modeling the technical system, the measurement data of 
the steady-state behavior of the system are of interest. 
10021] One refinement comprises reducing the number of 

measured values in a class by determining at least one 
representative value for the measured values of the class. 
Such a representative value may be: 

a) a mean value (for example a sliding mean value) of the 
measured values of the class, 

b) a maximum value of the measured values of the class, 

c) a minimum value of the measured values of the class, 

d) a median. 

[0028] In the case of variant d) , one advantage is that a 

value which actually exists can always be determined, whereas 
the mean value a) does not itself occur as a value. 

[0029] Depending on the application, a suitable choice can 

be made for determining the representative value of a class. 

[0030] An entire class with measured values may be picked 
out if it contains less than a predetermined number of 
measured values. 

[0031] Another refinement comprises that those measured 
values which vary by more than a predetermined threshold value 
from a predeterminable value are picked out . The 
predeterminable value may be a mean value of all the measured 
values of the class or a measured value to be expected in 
response to the respective setting parameters of the technical 
system. 

[0032] Within another development, the data obtained by 
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[means] way of designing are used for controlling a technical 
plant. In addition, the controlling of the technical plant 
can take place at the running time of the system-b — that is 
online . 

^ (i.e. , online) . 

[0033] Also specified for achieving the object is an 
arrangement for designing a technical system which has a 
processor unit [, which procoaaor unit io] set up [in] such -iet 
waY+ that measurement data of a predetermined system can be 
described on the basis of a substitute model. A numerical 
value for the quality of the substitute model can be 
determined by comparing the measurement data of the 
predetermined system with data determined by the substitute 
model. On the basis of the numerical value for the quality, 
the substitute model is adapted to be of as high a quality as 
possible. The substitute model adapted with regard to its 
quality is used for designing the technical system. 

[0034] This arrangement is particularly suitable for 

carrying out the method according to the invention or one of 
its developments explained above . 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0035] Exemplary embodiments of the invention are presented 
and explained below with reference to the drawing, 
9-in which: 

[figure 1 ahows] [0 03 6] Figure 1 is a block diagram which 
contains steps of a method for designing a technical system; 

[figure 2 ahov/g] [0037] Figure 2 is a schematic diagram of a 
recovery boiler; 

[figureo 3 5 show input variables, manipulated variables and 
output] [0038] Figure 3 is a chart showing the input v ariables 
of the recovery boiler-hr^-^. 
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[RcprcocntGd in figure 1 io] [0039] Figure 4 is a chart showing 
the manipulated variables of the recovery boiler; and 

[0040] Figure 5 is a chart showing the output variables of 

the recovery boiler. 

DETAILED DESCRIPTION OF THE INVENTION 
[0041] Figure 1 shows a block diagram which contains steps 

of a method for designing a technical system. In a step 101, 
a substitute model (preferably a regression model) is formed 
on the basis of measurement data. [Thia oubotitutc modol ia 
preferably a rcgrcaoion model . ] To adapt the substitute model 
created in step 101 to the measurement data [, that is] (i.e., 
to perform a refinement of the substitute model-hr+l so that 
the measurement data describe the substitute model in adequate 
approximation, a numerical value for the quality of the 
substitute model is determined in a step 102. This numerical 
value is determined by comparing the measurement data of the 
predetermined system with data determined by the substitute 
model . Each item of measurement data preferably receives a 
numerical value for the quality, which numerical value 
characterizes the deviation of the item of measurement data 
from the associated value determined by the substitute model. 
The sum of all the numerical values for the quality for all 
the measurement data determines an overall quality for the 
substitute model. In a step 103, the quality is maximized by 
minimizing the numerical value for the quality (or a negative 
quality for the coincidence of the substitute model with the 
predetermined system) . Once an appropriately high quality for 
the substitute model has been determined, this substitute 
model is used for designing the technical system in a step 
104. The designing may constitute both a new design (cf. step 
105) or an adaptation of an already existing technical system 
(cf . step 106) . 
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Figure 2 shows a schematic diagram of a recovery 
boiler. An exemplary embodiment of the method described above 
is illustrated below on the basis of the example of a 
"recovery boiler" . 

1004:3] In the paper and pulp industry, various chemicals 

and also heat and electrical energy are required for the 
digestion of pulp. The recovery boiler can be used to recover 
the chemicals used and also thermal energy from a concentrated 
process liquor (black liquor) . It is decisively 

significant to measure -te#9- the recovery of the chemicals 4ia 
of dcciaivc oignif icancc] for the cost-effectiveness of the 
overall plant. 

lOO^^l The black liquor is combusted in a char bed 201. As 

this happens, an alkali fusion is formed, flowing away via a 
line 202. In further process steps, the chemicals used are 
recovered from the constituents of the alkali fusion. 
Released heat of combustion is used for generating steam. The 
combustion of the waste liquor and consequently the recovery 
of the chemicals begins with the atomization of the black 
liquor via atomizer nozzles 204 into a combustion chamber 203. 
Particles of the atomized black liquor are dried as they fall 
through the hot flue gas. The dried liquor particles fall 
onto the char bed 201, with a first combustion and a chemical 
reduction taking place. Volatile constituents and reaction 
products pass into an oxidation zone, in which oxidizing 
reactions proceed and in which the combustion is completed. 
[Q0^5] Important targets for [the] controlling [of] the 

recovery boiler are the steam production for obtaining energy, 
the maintaining of emission values from environmental aspects 
and the efficiency of the chemical reduction. 
t0046] The combustion operation, and consequently the 
targets, are controlled in particular by the air supply on 
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three levels (Primary Air (PA) , Secondary Air (SA) , Tertiary 
Air (TA) ) . The overall process is subject to numerous 
influences, which have to be taken into account in the 
modeling : 

a) the measurement of the variables is often subject to 
strong f luctuat ions ; 

b) unmeasured and unmeasurable influencing variables 
exist ; 

c) any alteration of the parameters which can be set 
causes transient phenomena ; and 

d) the technical plant becomes soiled and is cleaned at 
predetermined intervals, [having the rcault] resulting 
each time [of] in a temporal drift in the system 
behavior . 

[0047] The measured variables of the overall process are 

divided into input variables (cf . [figure] Figure 3) and 
output variables (cf . [figure] Figure 5) . Measured values are 
stored every minute. Four of the input variables are at the 
same time also manipulated variables (also: parameters which 
can be set; cf . [figure] Figure 4) . The manipulated variables 
are [to be] regarded essentially as free parameters of the 
overall process which can be set independently of one another. 
Some of the other input variables are more or less dependent 
on the manipulated variables. According to one target, in the 
case of the recovery boiler_^ the "BL Front Pressure" and "BL 
Back Pressure" variables are always [to be] controlled such 
that they are the same. The four manipulated variables (cf. 
[figure] Figure 4) are preferably [to be] stored as 
manipulated variables (with the desired, preset value) and as 
input variables (with the measured, actual value) . 

t0048] In the case of the recovery boiler, one problem is 
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that of meeting targets which ; 1) are determined in dependence 
on the parameters to be set_j_ and 2)_ are defined by [moano of] 
measured variables. A three-stage procedure is chosen here to 
solve the problem: 

1. The targets to be considered are modeled by stochastic 
methods, these models being updated by new measurements 
(data-driven, empirical modeling) . For this step it is 
advisable to use not just a single model but instead global 
models for the identification of areas of interest in a 
parameter space determined by the targets and local models 
for the exact calculation of optimum operating points. The 
models used are assessed by measures of quality. 

2. If the models considered are not sufficiently accurate 
(measure of quality) [on account of] due to the data 
situation, new operating points are evaluated on a specific 
basis to improve the model (experimental design) . 
Furthermore, by the use of global stochastic optimizing 
methods with respect to the targets, attractive areas are 
identified in dependence on the current global model. 

3. For the local optimization, local models are devised and, if 
appropriate, the available sets of data are extended on a 
specific basis (experimental design) . 

[0049] The targets constitute physical -technical or 

business -management criteria, which generally have to conform 
to boundary conditions and/or safety conditions. Often a 
number of these criteria have to be considered simultaneously. 
The use of a stochastic model can be used in particular to 
simulate in the computer the target variables to be optimized 
and their dependence on the parameters to be set. This is 
necessary whenever measurements are very cost -intensive or 
very time-consuming. In the case of safety requirements, 
possible hazardous situations can be avoided. 
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[0050] In the case of the recovery boiler, an online 
optimization which is based on a number of data is necessary 
because the physical -chemical processes cannot be 
quantitatively modeled with sufficient accuracy and because 
the behavior of the plant is subject to fluctuations in the 
course of operation. Knowledge of this behavior must be 
continually extended by selective choice of new operating 
points. Therefore, the already described three-stage 
procedure of stochastic modeling and mathematical optimization 
is recommendable as part of online optimization. 

DESCRIPTION OF THE INPUT VARIABLES 
[0051] The a input variables (a e N, N: set of natural 

numbers) are generally dependent on n actuated variables n g N 
and on random effects. They can be described as follows: 

[0052] Let (Q,S, P) be a probability space and B"" be a 
Borel a algebra over 0"^ (D: set of real numbers) for each v e 
N. The input variables are represented by [means of] a x S 
- B^-measurable mapping cp : 

(p : X n -» 

[0053] The domain of the mapping (p is a Cartesian product 
of two sets. If the respective projections onto the 
individual sets are considered, the following mappings are 
obtained : 

[0054] 



(Px: ^ (o cp(x, co) for all X G 

[0055] 



(2) , 
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{p®:R^ -» R^, X -> <p{x, co) forall <o e Q (3) . 

[0056] [0057] is a stochastic 

{ process with an index set 

[0059] , and. a mapping cp® is, for each event co € Q, a 

path of this stochastic process. 

[0060] In the case of the recovery boiler, n=4 and a=14 

(after elimination of the "BL Back Pressure" variable) . 
[0061] On account of the required measurability of the 
mapping cpx, for each x € D"^ the mapping 9x is a random 
variable. Under suitable additional preconditions, 
expectation values and higher moments can be considered. This 
approach permits the step from stochastic models to 
deterministic optimization problems. In the case of a 
deterministic optimization problem, the target function can be 
set directly by [means of] a variable, whereas the stochastic 
variable influences the target function but does not permit a 
specific setting. 

DESCRIPTION OF THE OUTPUT VARIABLES 
[0062] The process model M of the recovery boiler is 

described as a function in dependence on the input variables 
and further random effects. In this case let (Q,S, P) be the 
above probability space. The process model M is then a x S 
- B^-measurable mapping: 
[0063] 



[0064] where b denotes the number of output variables. 
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[0065] Since the recovery boiler is subject to a cyclical 
temporal drift (from cleaning phase to cleaning phase) , a 
description with a time parameter is also conceivable. The 
output variables can be represented by B"" x S - B^-measurable 
mappings \\J : 
[0066] 

(x, co) ^ m((p(x, co), ©) (6) . 

[0067] If the respective projections onto the individual 

sets of the domain are considered, the following mappings are 
obtained 

[0068] 

y^f^iQ (0 ij/(x, £d) for all x e (7), 

vj/^iR'^ R^, X -> \\f{x,(o) for all w e Q (8). 

[0069] [0070] is a stochastic 

process with an index set 

[0072] D"", and the mapping \|;*' is, for each co e a path of 
this stochastic process. 

[0073] In the case of the recovery boiler, b=15 . 

[0074] The fact that in the definition of i|/ a distinction 

is not drawn between the events (o used does not mean that 
there is any restriction, since Q can be represented as a 
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Cartesian product of an Qi and an Q2 • The above 
representation consequently also comprises the model: 
[0075] 

(x, 0)1, C02) -> m((p(x, coi), 032) 



(9) , 

(10) . 



DESCRIPTION OF THE AVAILABLE SETS OF DATA 
[0076] With the descriptions in the two foregoing sections, 
the input variables and the output variables can be combined 
together to form measured variables (= measurement data) O. O 
is a X 5 - B™-measurable mapping with m = a + b and 
[0077] 



<U : R"" X R"^ (11) ^ 

Vv|/(x, cd). 



(^'^) ^ L../.. (12) 



[0078] If the respective projections onto the individual 

sets of the domain are again considered, the following 
mappings are obtained: 

Ox: Q R^, CO -> o(x, ©) for all x e R^ {13), 

<D^:R'' -> R^, X ^ O(x,(o) for all co e Q (14). 

-{- ^[0079] [0080] is a stochastic process 

with an index set 
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[0082] D"^, and the mapping is, for each © e Q, a path of 

this stochastic process. 

[0083] For each chosen related variable tupel x, in the 
case of the recovery boiler many realizations of Ox are 
determined and stored, i.e._^ for each Xj g [p- numerous 
realizations 
[0084] 

Ojk := cD(xj, (Djk) (15) 
with <o j k E Q ; k = 1,2, . . . , Vj ; 

V j e N ; j = 1,2, . . . , u ; u e N 

[0085] are considered. The stored sets of data Djk of the 
recovery boiler are consequently (n+m) tupels: 

[0086] 

^j^ = (oj^J ' ^ = 1-2 ; j = 1.2, ...,u (16) . 

[0087] In this case, Djiki is stored before Dj2k2 if 

Ol < 32) V (Oi = 32) A (ki < k2)) 
[0088] ^ 

[0089] applies. 

DATA COMPRESSION BY CLASSIFICATION OF THE PARAMETERS 
[0090] Since for each manipulated variable tupel x there 
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are generally a number of realizations of Ox, a classification 
of the parameters by forming arithmetic mean values is a 
suitable operation as a first step of the statistical data 
analysis on account of the complex stochastic properties of 
the process to be considered. Moreover, sets of data which 
are obviously erroneous are picked out. A set of data which 
is obviously erroneous is, for example, a physically 
impossible measurement which cannot occur in reality, -tiii 
particular on account of] particularly due to a setting made. 

[0091] Procedure : 

1. Sets of data in which the "BL Front Pressure" variable is 
not equal to the "BL Back Pressure" variable are picked out, 
since these two values must be equal according to the 
default of the plant control. The data loss is very small. 

2. The sets of data are divided into classes in which the four 
setting parameters (PA, SA, TA, BL Front Pressure, see 
above) are constant in temporal succession, 

[0092] i.e.^ the j th class comprises the data sets Dj, . 

3. Classes in which there are fewer than 3 0 data sets are 
picked out in order that transient phenomena have no great 
influence . 

4. For each class, an arithmetic mean value <l)-|~and an empirical 
standard deviation Sj are determined for all the measured 
variables : 

[0093] 
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f 




V k = l 



(18) . 



1 




5. Classes in which the mean values for the variables PA, SA, 
TA or BL Front Pressure are too far away from the 
corresponding setting parameters are picked out. Therefore, 
in these classes it [was] is not possible to reach the 
setting values. 

STATISTICAL CHARACTERISTIC VARIABLES FOR THE GIVEN 
CLASSES AND THEIR GRAPHIC REPRESENTATION 

[0094] Apart from the arithmetic mean values and the 
empirical standard deviations which were determined for the 
individual classes, a common standard deviation s is also 
determined according to 

[0095] 
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1^ 



j=l_ 



j=i 



(19) 



[0096] where u stands for the number of classes (here 205) 

and V stands for the sum of Vj , i.e._j_ v is the number of all 
the measured values used (here 38915) . 

LINEAR REGRESSION MODELS FOR FUNCTION APPROXIMATIONS 
[0097] For each measured variable (item of measurement 

data) O^^^ (1=1,2, . . . ,m) , a linear regression model is 
calculated on the basis of the arithmetic mean over the 
classes in dependence on the quadratic combination of the four 
setting parameters. In the following representation, x e R^, 
where 

[0098] x^^' : Primary Air (PA) 

[0099] x^^' : Secondary Air (SA) 

[00100] x^^^ : Tertiary Air (TA) 

[00101] x^^^ : Black Liquor (BL) Front Pressure 

[00102] applies. u e N denotes the number of classes. Each 
measured variable O*^^ is modeled by 
[00103] 



aj^r(x) + eiCco) 



(20) 



[00104] with ai e R^^ . Here the following applies: 
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[00105] 



15 



^1^2. ^1^3' ^1^4. «4)' 



(21) 



(22) , 



[00106] i.e._^ polynomials of the second degree are adapted 
to the measurement data, and 

[00107] 



ei : Q ^ R 

[00108] is a random variable with the expected value 0. 
[00109] The vector ai is determined by the method of 
[00110] least squares, but the -h ^[00111] 



(23) 



arithmetic means 



[00113] are used instead of the 
original data sets 



[00116] This procedure is suitable, since linear regression 
models estimate in particular expected values. This results 
in the following minimization problem: 

[00117] 
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(24) . 



[00118] Let ai be the optimum point of the quadratic 
minimization problem from equation (24) . Furthermore, the 
following applies: 

[00119] 













^(15) 



(25) , 



1 ^ _(. 



(26) 



j = l 



[00120] To validate the regression theorem, a coefficient of 
determination is calculated according to 

[with 



^ [00121] 



-T- —2 



(27) 



[00122] 
[00123] 
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[00124] The closer R^i is to 1, the better the dependent 
variable is represented by the independent variables 
(o < < l)| 

[00125] In addition, a maximum E^^-'^ax for an absolute value 
of the deviation of the data from the model is specified as 
[00126] 



max 



- yl 



(29) . 



[00127] is that value below which at least 90% of the 

absolute values of the deviations of the data from the model 
lie. By analogy with this, E go% is that value below which at 
least 80% of the absolute values of the deviations of the data 
from the model lie. With the optimum point ai of the 
minimization problem according to equation (24) , a model O'^^ 
of the expected value~of the measured variable O'"-^ can be 
specified as 

[00128] 



O^^^ := R (30), 



X af r(x) 



(31) . 



[00129] In particular, the gradient VO^^^ J can be 
analytically specified by 
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[00130] 



Vd^%) = ;^(x) • ai for all x € (32 

[00131] The above-described method and system are 

illustrative of the principles of the present invention. 
Nxjmerous modifications and adaptations will be readily- 
apparent to those skilled in this art without departing from 
the spirit and scope of the present invention. 
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ABSTRACT 

[Method and arrangcTnont for doDigning a tGchnical oyatcm 
To achieve the object, a method ] [00132] A method is 

provided for designing a technical system in which measurement 
data of a predetermined system are described based on [the 
baoio of] a substitute model [is □pccificd] . A numerical 
value for the quality of the substitute model is determined by 
comparing the measurement data of the predetermined system 
with data determined by the substitute model. On the basis of 
the numerical value for the quality, the substitute model is 
adapted to be of as high a quality as possible. The 
substitute model adapted with regard to its quality is used 
for designing the technical system. 
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SPECIFICATION 
TITLE 

METHOD AND ARRANGEMENT FOR DESIGNING A TECHNICAL SYSTEM 
BACKGROUND OF THE INVENTION 
Field of the Invention 

[0001] The invention relates to a method and arrangement 

for designing a technical system. 

Description of the Related Art 

[0002] The system behavior of a technical system, for 
example, a process engineering plant or system in heavy 
industry, depends on numerous parameters. In the course of 
designing such a system, particularly in the case of a new 
design or when adapting or adjusting an already existing 
system, it is necessary to comply with preconditions such as 
the cost-ef f ectiveness or environmental impact of the system. 
Each precondition is formulated as a target function, the 
general aim being to optimize this with respect to the other 
target functions . 

SUMMARY OF THE INVENTION 
[00 03] The object of the invention is to make it possible 

for a technical system to be designed on the basis of 
measurement data of a predetermined system. Specifically with 
regard to the optimization of the existing system or with 
regard to an optimized new design of a system, such use of 
known measurement data is of great significance. 

[0004] This object is achieved a method and apparatus 
described below. 

[00 05] To achieve the object, the present invention 
provides a method for designing a technical system in which 
measurement data of a predetermined system are described on 
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the basis of a substitute model. A numerical value for the 
quality of the substitute model is determined by comparing the 
measurement data of the predetermined system with data 
determined by the substitute model. Based on the numerical 
value for the quality, the substitute model is adapted to be 
of as high a quality as possible. 

[0006] The substitute model adapted with regard to its 
quality is used for designing the technical system. 

[0007] The measurement data obtained from many different 
realized systems are used for describing the substitute model. 
The substitute model attempts to replicate the predetermined 
system as well as possible. The numerical value for the 
quality of the replication is determined by comparing the 
actual measurement data with the data obtained on the basis of 
the substitute model. A great difference between the 
measurement data and the data of the substitute model 
corresponds to poor quality, that is, a poor mapping of the 
predetermined system into the substitute model. The numerical 
value for the quality is used to adapt the substitute model to 
make the quality itself become as high as possible and 
consequently to make the substitute model describe the 
predetermined system as well as possible. The high-quality 
substitute model obtained in this way is used for designing 
the technical system. 

[0008] Designing is understood in a general sense as 
meaning both the new design of a technical system and the 
adaptation of an already existing technical system. 

[0009] One development comprises that the substitute model 
is a regression model. The regression model is based on the 
description 



[0010] 



Yi = fp(Xi) + ei 



[0011] 



where 



[0012] 



(Yi, Xi) 



denotes predetermined pairs of 
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values, 

[0013] fp denotes a function which is dependent on a 
parameter p, and 

[0014] ei denotes an error. 
[0015] It is then intended to minimize the error (as a 

function of P) : 

n 

Z4 = <p(p) . 

i = l 

[0016] 

[0017] If the following example 
[0018] Y = Po + Pix + p2x2 + e 

[0019] is taken as a basis, the functional relationship is 
of a quadratic order, while the regression model (function, 
dependent on P) is linear. 

[0020] In another development, the quality can be determined 

on the basis of a mean square deviation of the measurement data 
from the data determined by the siobstitute model . The adaptation 
of the substitute model takes place by minimizing the mean square 
deviation. 

[0021] One refinement comprises that the measurement data are 

sorted according to their quality, with respect to the deviation 
of the latter from the data determined by the substitute model, 
and a predetermined number of n% of the worst measurement data 
are picked out. Consequently, a quality is determined for each 
item of measurement data, the volume of measurement data, 
preferably in the form of a list, being sorted according to their 
quality and the n% of the worst measurement data, or the n worst 
measurement data, being picked out. In particular, it must be 
checked whether the n% of worst measurement data or n worst 
measurement data lie in a continuous range. If this is the case, 
these measurement data are not picked out, since they vezy 
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probably do not represent measurement errors but a continuous 
range which is not mapped accurately enough by the substitute 
model . 

[0022] Another development comprises that the measurement 
data are subjected to a preprocessing operation. Since, in an 
actual predetermined system, a large amount of measurement 
data occur per unit of time, it is advisable to subject these 
measurement data to a preprocessing operation and consequently 
ensure that largely significant measurement data are taken for 
forming the substitute model. The preprocessing operation 
preferably finds its expression in a reduction in the number 
of measurement data. 

[0023] In the process, the measurement data are divided 
into classes according to predetermined criteria. The 
measured values of a class are assessed and those measured 
values for which the assessment lies below a predetermined 
first threshold value are picked out. The picking out of the 
measured values results in a reduction with regard to the 
number of measured values. This results in a significantly 
reduced number of measured values for a further processing 
operation. The further processing operation can take place 
with less computing effort in comparison with the unreduced 
number of measured values . 

[0024] The classes themselves can also be assessed. In 
particular, a class for which the assessment lies below a 
predetermined second threshold value can be picked out . As a 
result, an additional reduction of the number of measured 
values is obtained. 

[0025] Another development of the preprocessing operation 
comprises that a criterion for the classification is that, for 
each class, measured values are determined as a default for 
setting parameters of the technical system. The technical 
system is typically set on the basis of a predetermined number 
of setting parameters; after setting, a reaction of the system 

- 4 - Substitute Specification 



to the setting parameters takes place (usually with a time 
delay) (transient response, transient phenomenon of the 
system) . After setting, consequently a certain set of 
measured values which can be assigned to the transient 
phenomenon are recorded, with measured values continuing to 
occur after the transient phenomenon has come to an end 
(transition to steady-state operation) , and are assigned to 
the predetermined set of setting parameters. Adjusting the 
setting parameters determines a new class. All of the 
measured values which respectively occur after adjustment of 
the setting parameters belong in a class of their own. 

[0026] In addition, measured values of a class which can be 
assigned to the respective transient phenomenon can be picked 
out. Furthermore, erroneous measured values can be picked 
out. The setting of large technical systems is in many cases 
directed at long-term steady-state operation. It is advisable 
for the measured values which relate to the transient 
phenomenon (of short duration in relation to the steady- state 
operation after the transient phenomenon has come to an end) 
to be picked out, since they have the effect of falsifying 
measured values for steady-state operation. In particular, 
when modeling the technical system, the measurement data of 
the steady-state behavior of the system are of interest. 

[0027] One refinement comprises reducing the number of 
measured values in a class by determining at least one 
representative value for the measured values of the class. 
Such a representative value may be: 

a) a mean value (for example a sliding mean value) of the 
measured values of the class, 

b) a maximum value of the measured values of the class, 

c) a minimum value of the measured values of the class, 

d) a median. 

[0028] In the case of variant d) , one advantage is that a 
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value which actually exists can always be determined, whereas 
the mean value a) does not itself occur as a value. 

[0029] Depending on the application, a suitable choice can 

be made for determining the representative value of a class. 

[0030] An entire class with measured values may be picked 
out if it contains less than a predetermined number of 
measured values. 

[0031] Another refinement comprises that those measured 
values which vary by more than a predetermined threshold value 
from a predeterminable value are picked out . The 
predeterminable value may be a mean value of all the measured 
values of the class or a measured value to be expected in 
response to the respective setting parameters of the technical 
system. 

[0032] Within another development, the data obtained by way 
of designing are used for controlling a technical plant. In 
addition, the controlling of the technical plant can take 
place at the running time of the system (i.e., online). 

[0033] Also specified for achieving the object is an 
arrangement for designing a technical system which has a 
processor unit set up such that measurement data of a 
predetermined system can be described on the basis of a 
substitute model. A numerical value for the quality of the 
substitute model can be determined by comparing the 
measurement data of the predetermined system with data 
determined by the substitute model. On the basis of the 
numerical value for the quality, the substitute model is 
adapted to be of as high a quality as possible. The 
substitute model adapted with regard to its quality is used 
for designing the technical system. 

[0034] This arrangement is particularly suitable for 
carrying out the method according to the invention or one of 
its developments explained above. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[0035] Exemplary embodiments of the invention are presented 
and explained below with reference to the drawing, in which: 
[0036] Figure 1 is a block diagram which contains steps of 
a method for designing a technical system; 

[003 7] Figure 2 is a schematic diagram of a recovery 

boiler; 

[0038] Figure 3 is a chart showing the input variables of 

the recovery boiler; 

[0 03 9] Figure 4 is a chart showing the manipulated 
variables of the recovery boiler; and 

[0040] Figure 5 is a chart showing the output variables of 
the recovery boiler. 

DETAILED DESCRIPTION OF THE INVENTION 
[0041] Figure 1 shows a block diagram which contains steps 
of a method for designing a technical system. In a step 101, 
a substitute model (preferably a regression model) is formed 
on the basis of measurement data. To adapt the substitute 
model created in step 101 to the measurement data (i.e., to 
perform a refinement of the substitute model) so that the 
measurement data describe the substitute model in adequate 
approximation, a numerical value for the quality of the 
substitute model is determined in a step 102 . This numerical 
value is determined by comparing the measurement data of the 
predetermined system with data determined by the substitute 
model. Each item of measurement data preferably receives a 
numerical value for the quality, which numerical value 
characterizes the deviation of the item of measurement data 
from the associated value determined by the substitute model. 
The sum of all the numerical values for the quality for all 
the measurement data determines an overall quality for the 
substitute model. In a step 103, the quality is maximized by 
minimizing the numerical value for the quality (or a negative 
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quality for the coincidence of the substitute model with the 
predetermined system) . Once an appropriately high quality for 
the substitute model has been determined, this substitute 
model is used for designing the technical system in a step 
104. The designing may constitute both a new design (cf. step 
105) or an adaptation of an already existing technical system 
(cf . step 106) . 

[0 042] Figure 2 shows a schematic diagram of a recovery 

boiler. An exemplary embodiment of the method described above 
is illustrated below on the basis of the example of a 
"recovery boiler". 

[0043] In the paper and pulp industry, various chemicals 

and also heat and electrical energy are required for the 
digestion of pulp. The recovery boiler can be used to recover 
the chemicals used and also thermal energy from a concentrated 
process liquor (black liquor) . It is decisively significant 
to measure the recovery of the chemicals for the cost- 
effectiveness of the overall plant. 

[0044] The black liquor is combusted in a char bed 201. As 
this happens, an alkali fusion is formed, flowing away via a 
line 202. In further process steps, the chemicals used are 
recovered from the constituents of the alkali fusion. 
Released heat of combustion is used for generating steam. The 
combustion of the waste liquor and consequently the recovery 
of the chemicals begins with the atomization of the black 
liquor via atomizer nozzles 204 into a combustion chamber 203. 
Particles of the atomized black liquor are dried as they fall 
through the hot flue gas. The dried liquor particles fall 
onto the char bed 201, with a first combustion and a chemical 
reduction taking place. Volatile constituents and reaction 
products pass into an oxidation zone, in which oxidizing 
reactions proceed and in which the combustion is completed. 

[0045] Important targets for controlling the recovery 
boiler are the steam production for obtaining energy, the 
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maintaining of emission values from environmental aspects and 
the efficiency of the chemical reduction, 

[0046] The combustion operation, and consequently the 

targets, are controlled in particular by the air supply on 
three levels (Primary Air (PA) , Secondary Air (SA) , Tertiary 
Air (TA) ) . The overall process is subject to numerous 
influences, which have to be taken into account in the 
modeling : 

a) the measurement of the variables is often subject to 
strong fluctuations; 

b) unmeasured and unmeasurable influencing variables 
exist ; 

c) any alteration of the parameters which can be set 
causes transient phenomena; and 

d) the technical plant becomes soiled and is cleaned at 
predetermined intervals, resulting each time in a 
temporal drift in the system behavior. 

[0047] The measured variables of the overall process are 
divided into input variables (cf. Figure 3) and output 
variables (cf. Figure 5). Measured values are stored every 
minute. Four of the input variables are at the same time also 
manipulated variables (also: parameters which can be set; cf . 
Figure 4) . The manipulated variables are regarded essentially 
as free parameters of the overall process which can be set 
independently of one another. Some of the other input 
variables are more or less dependent on the manipulated 
variables. According to one target, in the case of the 
recovery boiler, the "BL Front Pressure" and "BL Back 
Pressure" variables are always controlled such that they are 
the same. The four manipulated variables (cf . Figure 4) are 
preferably stored as manipulated variables (with the desired, 
preset value) and as input variables (with the measured, 
actual value) . 
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[0048] In the case of the recovery boiler, one problem is 

that of meeting targets which: 1) are determined in dependence 
on the parameters to be set, and 2) are defined by measured 
variables. A three-stage procedure is chosen here to solve 
the problem: 

1. The targets to be considered are modeled by stochastic 
methods, these models being updated by new measurements 
{data-driven, empirical modeling) . For this step it is 
advisable to use not just a single model but instead global 
models for the identification of areas of interest in a 
parameter space determined by the targets and local models 
for the exact calculation of optimum operating points. The 
models used are assessed by measures of quality. 

2. If the models considered are not sufficiently accurate 
(measure of quality) due to the data situation, new 
operating points are evaluated on a specific basis to 
improve the model (experimental design) . Furthermore, by 
the use of global stochastic optimizing methods with respect 
to the targets, attractive areas are identified in 
dependence on the current global model . 

3 . For the local optimization, local models are devised and, if 
appropriate, the available sets of data are extended on a 
specific basis (experimental design) . 

[0049] The targets constitute physical -technical or 

business -management criteria, which generally have to conform 
to boundary conditions and/or safety conditions. Often a 
number of these criteria have to be considered simultaneously. 
The use of a stochastic model can be used in particular to 
simulate in the computer the target variables to be optimized 
and their dependence on the parameters to be set. This is 
necessary whenever measurements are very cost -intensive or 
very time-consuming. In the case of safety requirements, 
possible hazardous situations can be avoided. 

[0050] In the case of the recovery boiler, an online 
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optimization which is based on a number of data is necessary 
because the physical -chemical processes cannot be 
quantitatively modeled with sufficient accuracy and because 
the behavior of the plant is subject to fluctuations in the 
course of operation. Knowledge of this behavior must be 
continually extended by selective choice of new operating 
points. Therefore, the already described three- stage 
procedure of stochastic modeling and mathematical optimization 
is recommendable as part of online optimization. 

DESCRIPTION OF THE INPUT VARIABLES 

[0051] The a input variables (a e N, N: set of natural 
numbers) are generally dependent on n actuated variables n e N 
and on random effects. They can be described as follows: 

[0052] Let (Q,S, P) be a probability space and B'' be a 
Borel <j algebra over D"" {□: set of real numbers) for each v € 
N. The input variables are represented by a B"^ x S - B*- 
measurable mapping <p : 



cp : R" X Q 



(1) 



[0053] The domain of the mapping cp is a Cartesian product 
of two sets. If the respective projections onto the 
individual sets are considered, the following mappings are 
obtained: 

(Px:^ R^f ca <p(x, co) for all x e (2) 

[0054] 

[0055] R^ X <p(x,co) for all oj e Q (3) 

f [0057] is a stochastic process 

[0 056] ^ with an index set 

[0058] U^, and a mapping 9^ is, for each event 00 e Q, a 

path of this stochastic process. 
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[0059] In the case of the recovery boiler, n=4 and a=14 

(after elimination of the "BL Back Pressure" variable) . 
[OOSO] On account of the required measurability of the 
mapping cpx, for each x s D"^ the mapping (px is a random 
variable. Under suitable additional preconditions, 
expectation values and higher moments can be considered. This 
approach permits the step from stochastic models to 
deterministic optimization problems. In the case of a 
deterministic optimization problem, the target function can be 
set directly by a variable, whereas the stochastic variable 
influences the target function but does not permit a specific 
setting . 

DESCRIPTION OF THE OUTPUT VARIABLES 

[0061] The process model M of the recovery boiler is 
described as a function in dependence on the input variables 
and further random effects. In this case let (Q,S, P) be the 
above probability space. The process model M is then a x S 
- B"^ -measurable mapping: 

M : X Q -> ^A^ 

[0062] 

[0 063] where b denotes the number of output variables. 

[0064] Since the recovery boiler is subject to a cyclical 
temporal drift {from cleaning phase to cleaning phase) , a 
description with a time parameter is also conceivable. The 
output variables can be represented by B'^ x S - B^-measurable 
mappings \|/ : 

V : R" X Q R^ (5), 

[0065] ^""'"^^ "W^' (6) . 

[0066] If the respective projections onto the individual 

sets of the domain are considered, the following mappings are 
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obtained 



[0067] 



\j/x: n R , G> co) for all x e r" 

M/®:R" -> R^, X vk(x, ©) for all co e Q 

[0069] is a stochastic process 



(7) , 



with an index set 



[0068] 

[0 070] and the mapping ii/*" is, for each co e Q, a path of 

this stochastic process. 

[0071] In the case of the recovery boiler, b=15. 

[0072] The fact that in the definition of y a distinction 
is not drawn between the events (o used does not mean that 
there is any restriction, since Q can be represented as a 
Cartesian product of an Clx and an Q2 . The above 
representation consequently also comprises the model : 

V : R^ X ni X Q2 ( 9 ) , 

(x, Q>i, 102) -> m(cp(x, coi), ©2) (10) 

[0073] 

DESCRIPTION OF THE AVAILABLE SETS OF DATA 

[0074] With the descriptions in the two foregoing sections, 
the input variables and the output variables can be combined 
together to form measured variables (= measurement data) O. 
<I> is a B° X S - B'^-measurable mapping with m = a + b and 

O : R" X Q R"' (11) 

. V f <p(x, a>)) 

[0075] 

[0076] If the respective projections onto the individual 
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sets of the domain are again considered, the following 
mappings are obtained: 



Ox: ^2 -> R^, 0, <d(x, o>) for all x e 



(13) , 



0<»:R" R"", X <I>(x,a)) for all o> € Q (14) . 

( [0078] is a stochastic process 

Ov/ X € 

[0077] with an index set 

[0079] and the mapping O"" is, for each a e Q, a path of 

this stochastic process. 

[0080] For each chosen related variable tupel x, in the 
case of the recovery boiler many realizations of Ox are 
determined and stored, i.e., for each Xj e numerous 
realizations 

with o) j ^. e n ; k = 1,2, . . . , Vj ; 
Vj € N ; j = 1,2, . . . , u ; u € N 

[0081] 

[0082] are considered. The stored sets of data Djk of the 

recovery boiler are consequently (n+m) tupels: 



- (^''^J. k = 1,2 Vj; j = 1,2,, 

[0083] 

[0084] In this case, Djiki is stored before Dj2k2 if 

[0085] ^ " = ^ ^2)) 

[0086] applies. 

DATA COMPRESSION BY CLASSIFICATION OF THE PARAMETERS 
[0087] Since for each manipulated variable tupel x there 
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are generally a number of realizations of Ox, a classification 
of the parameters by forming arithmetic mean values is a 
suitable operation as a first step of the statistical data 
analysis on account of the complex stochastic properties of 
the process to be considered. Moreover, sets of data which 
are obviously erroneous are picked out. A set of data which 
is obviously erroneous is, for example, a physically 
impossible measurement which cannot occur in reality, 
particularly due to a setting made. 

[0088] Procedure: 

1. Sets of data in which the "BL Front Pressure" variable is 
not equal to the "BL Back Pressure" variable are picked out, 
since these two values must be equal according to the 
default of the plant control. The data loss is very small. 

2. The sets of data are divided into classes in which the four 
setting parameters (PA, SA, TA, BL Front Pressure, see 
above) are constant in temporal succession, 

[0089] i.e., the j th class comprises the data sets Dj,. 

3. Classes in which there are fewer than 3 0 data sets are 
picked out in order that transient phenomena have no great 
influence . 

4. For each class, an arithmetic mean value Oj^and an empirical 
standard deviation Sj are determined for all the measured 
variables : 
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(17) 




Sj = 



(18) 



1 



( 



2h 



1 



of-; 




- 1 



k=l 



[0090] 



5. Classes in which the mean values for the variables PA, SA, 
TA or BL Front Pressure are too far away from the 
corresponding setting parameters are picked out. Therefore, 
in these classes it is not possible to reach the setting 
values . 

STATISTICAL CHARACTERISTIC VARIABLES FOR THE GIVEN CLASSES AND 
THEIR GRAPHIC REPRESENTATION 

[0091] Apart from the arithmetic mean values and the 
empirical standard deviations which were determined for the 
individual classes, a common standard deviation s is also 
determined according to 



and V stands for the sum of Vj , i.e., v is the number of all 
the measured values used (here 38915) . 




(19) 



[0092] 



[0093] 



where u stands for the number of classes (here 2 05) 
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LINEAR REGRESSION MODELS FOR FUNCTION APPROXIMATIONS 
[0094] For each measured variable (item of measurement 
data) O^^^ (i=l , 2 , . . . , m) , a linear regression model is 
calculated on the basis of the arithmetic mean over the 
classes in dependence on the quadratic combination of the four 
setting parameters. In the following representation, x e R'*, 



where 








[0095] 




Primary Air (PA) 




[0096] 


x^^>: 


Secondary Air (SA) 




[0097] 


x«^>: 


Tertiary Air (TA) 




[0098] 


x«^»: 


Black Liquor (BL) Front Pressure 




[0099] 


applies . 


u e N denotes the number of classes. 


Each 


measured 


variable O 


is modeled by 




[00100] 


0^i)(x, o) ■ 


= ajr(x) + ei{Q>) 


(20) 


[00101] 


with ai e 


R"""^ . Here the following applies: 






r : 


r15 


(21) 



fa^ C4f ^ ( 1, ;i, ^2' ^3' ^4. 4r ^1 

;iC2.. ;i^3. <;i^4. ^2^3. ^2C4. ';3^4 ) (22 ) , 

[00102] 

[00103] i.e., polynomials of the second degree are adapted 
to the measurement data, and 

ei : Q R (23) 

[00104] 

[00105] is a random variable with the expected value 0. 
[0 0106] The vector ai is determined by the method of 

[00107] least squares, but the y- 

arithmetic means 1^1' )• 

[00108] ^ 
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[00109] are used instead of the 
original dat^ sets 



[00110] 



[00111] This procedure is suitable, since linear regression 
models estimate in particular expected values. This results 
in the following minimization problem: 



,15 



(24) . 



[00112] 

[00113] Let ai be the optimum point of the quadratic 
minimization problem from equation (24) . Furthermore, the 
following applies: 



(25) , 



Yi : = Jo^^i) e R (26) 

j = l 

To validate the regression theorem, a coefficient of 



[00114] 
[00115] 

determination is calculated according to 

\2 



[00116] 
[00117] 



2 ^(i)T^i) _ 



(27) 



[00118] 



^(i) 



(28) . 
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[00119] The closer R^i is to 1, the better the dependent 
variable is represented by the independent variables 
(o < R? < l) 

[00120] In addition, a maximum E'^^rnax for an absolute value 
of the deviation of the data from the model is specified as 

[00121] 

[00122] is that value below which at least 90% of the 

absolute values of the deviations of the data from the model 
lie. By analogy with this, E*'-'8o% is that value below which at 
least 80% of the absolute values of the deviations of the data 
from the model lie. With the optimum point ai of the 
minimization problem according to equation (24) , a model O^^' 
of the expected value~f the measured- variable O^^' can be 
specified as 



[00123] 



ai^r(x) (31) 



[00124] In particular, the gradient VO^^' can be analytically 
specified by 

va^%) = ^(x) • ai for ail x e (32 
dx 

[00125] 

[00126] The above -described method and system are 

illustrative of the principles of the present invention. 
Numerous modifications and adaptations will be readily 
apparent to those skilled in this art without departing from 
the spirit and scope of the present invention. 
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ABSTRACT 

[00127] A method is provided for designing a technical 
system in which measurement data of a predetermined system are 
described based on a substitute model . A numerical value for 
the quality of the substitute model is determined by comparing 
the measurement data of the predetermined system with data 
determined by the substitute model. On the basis of the 
numerical value for the quality, the substitute model is 
adapted to be of as high a quality as possible. The 
substitute model adapted with regard to its quality is used 
for designing the technical system. 
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Description 



Method and arrangement for designing a technical system 

The invention relates to a method and 
arrangement for designing a technical system. 

The system behavior of a technical system, for 
example a process engineering plant or system in heavy 
industry, depends on numerous parameters. In the 
course of designing such a system, that is in 
particular in the case of a new design or when adapting 
or adjusting an already existing system, it is 
necessary to comply with preconditions, for example 
with regard to the cost-effectiveness or environmental 
impact of the system. Each precondition is formulated 
as a target function, the optimization of which with 
regard to the other target functions is the general 
aim. 

The object of the invention is to make it 
possible for a technical system to be designed on the 
basis of measurement data of a predetermined system. 
Specifically with regard to the optimization of the 
existing system or with regard to an optimized new 
design of a system, such use of known measurement data 
is of great significance. 

This object is achieved according to the 
features of the independent patent claims. 
Developments of the invention also emerge from the 
dependent claims . 

To achieve the object, a method for designing a 
technical system in which measurement data of a 
predetermined system are described on the basis of a 
substitute model is specified. A numerical value for 
the quality of the substitute model is determined by 
comparing the measurement data of the predetermined 
system with data determined by the substitute model. 
On the basis of the numerical value for the quality, 
the substitute model is adapted to be of as high a 
quality as possible. 
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The substitute model adapted with regard to its quality 
is used for designing the technical system. 

The measurement data obtained from many 
different realized systems are used for describing the 
5 substitute model. With the substitute model it is 
attempted to replicate the predetermined system as well 
as possible. The numerical value for the quality of 
the replication is determined by comparing the actual 
measurement data with the data obtained on the basis of 

10 the substitute model. A great difference between the 
measurement data and the data of the substitute model 
corresponds to poor quality, that is poor mapping of 
the predetermined system into the substitute model. 
The numerical value for the quality is used to adapt 

15 the substitute model to make the quality itself become 
as high as possible and consequently to make the 
substitute model describe the predetermined system as 
well as possible. The high-quality substitute model 
obtained in this way is used for designing the 

20 technical system. 

Designing is understood in a general sense as 
meaning both the new design of a technical system and 
the adaptation of an already existing technical system. 

One development comprises that the substitute 

25 model is a regression model. The regression model is 
based on the description 

Yi = fp{Xi) + ei 

where 

(Yir Xi) denotes predetermined pairs of values, 
30 f(5 denotes a function which is dependent on 

a parameter p, and 
ei denotes an error. 
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It is then intended to minimize the error (as a 
function of P) : 

n 
i = l 

If the following example 
Y = po + Pix + Psx^ + e 

is taken as a basis, the functional relationship is of 
a quadratic order, while the regression model 
(function, dependent on p) is linear. 

In another development, the quality can be 
determined on the basis of a mean square deviation of the 
measurement data from the data determined by the 
substitute model. The adaptation of the substitute model 
takes place by minimizing the mean square deviation. 

One refinement comprises that the measurement data 
are sorted according to their quality, with respect to the 
deviation of the latter from the data determined by the 
substitute model, and a predetermined number of n% of the 
worst measurement data are picked out. Consequently, a 
quality is determined for each item of measurement data, 
the volume of measurement data, preferably in the form of 
a list, being sorted according to their quality and the n% 
of the worst measurement data, or the n worst measurement 
data, being picked out. In particular, it must be checked 
whether the n% of worst measurement data or n worst 
measurement data lie in a continuous range. If this is 
the case, these measurement data are not picked out, since 
they very probably do not represent measurement errors but 
a continuous range which is not mapped accurately enough 
by the substitute model. 

Another development comrises that the measure- 
ment data are subjected to a preprocessing operation. 
Since, in an actual predetermined system, a large 
amount of measurement data occur per 
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unit of time, it is advisable to subject these 
measurement data to a preprocessing operation and 
consequently ensure that largely significant 
measurement data are taken for forming the substitute 
5 model. The preprocessing operation preferably finds 
its expression in a reduction in the number of 
measurement data. 

In the process, the measurement data are 
divided into classes according to predetermined 

10 criteria. The measured values of a class are assessed 
and those measured values for which the assessment lies 
below a predetermined first threshold value are picked 
out. The picking out of the measured values results in 
a reduction with regard to the number of measured 

15 values. This results in a significantly reduced number 
of measured values for a further processing operation. 
The further processing operation can take place with 
less computing effort in comparison with the unreduced 
number of measured values . 

20 The classes themselves can also be assessed. 

In particular, a class for which the assessment lies 
below a predetermined second threshold value can be 
picked out. As a result, an additional reduction of 
the number of measured values is obtained. 

25 Another development of the preprocessing 

operation comprises that a criterion for the 
classification is that, for each class, measured values 
are determined as a default for setting parameters of 
the technical system. The technical system is 

30 typically set on the basis of a predetermined number of 
setting parameters; after setting, a reaction of the 
system to the setting parameters takes place (usually 
with a time delay) (transient response, transient 
phenomenon of the system) . After setting, consequently 

35 a certain set of measured values which can be assigned 
to the transient phenomenon are recorded, with measured 
values continuing to occur after the transient 
phenomenon has come to an end (transition to steady- 
state operation) , and 

40 
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are assigned to the predetermined set of setting 
parameters. By adjusting the setting parameters, a new 
class is determined. All the measured values which 
respectively occur after adjustment of the setting 
5 parameters belong in a class of their own. 

In addition, measured values of a class which 
can be assigned to the respective transient phenomenon 
can be picked out. Furthermore, erroneous measured 
values can be picked out. The setting of large 

10 technical systems is in many cases directed at long- 
term steady-state operation. It is advisable for the 
measured values which relate to the transient 
phenomenon (of short duration in relation to the 
steady-state operation after the transient phenomenon 

15 has come to an end) to be picked out, since they have 
the effect of falsifying measured values for steady- 
state operation. In particular when modeling the 
technical system, the measurement data of the steady- 
state behavior of the system are of interest. 

2 0 One refinement comprises reducing the number of 

measured values in a class by determining at least one 
representative value for the measured values of the 
class. Such a representative value may be: 

a) a mean value (for example a sliding mean value) 
25 of the measured values of the class, 

b) a maximum value of the measured values of the 
class , 

c) a minimum value of the measured values of the 
class , 

30 d) a median. 

In the case of variant d) , one advantage is 

that a value which actually exists can always be 

determined, whereas the mean value a) does not itself 

occur as a value. 
35 Depending on the application, a suitable choice 

can be made for determining the representative value of 

a class. 
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An entire class with measured values may be 
picked out if it contains less than a predetermined 
number of measured values. 

Another refinement comprises that those 
5 measured values which vary by more than a predetermined 
threshold value from a predeterminable value are picked 
out . The predeterminable value may be a mean value of 
all the measured values of the class or a measured 
value to be expected in response to the respective 
10 setting parameters of the technical system. 

Within another development, the data obtained 
by means of designing are used for controlling a 
technical plant. In addition, the controlling of the 
technical plant can take place at the running time of 
15 the system, that is online. 

Also specified for achieving the object is an 
arrangement for designing a technical system which has 
a processor unit, which processor unit is set up in 
such a way that measurement data of a predetermined 
20 system can be described on the basis of a substitute 
model. A numerical value for the quality of the 
substitute model can be determined by comparing the 
measurement data of the predetermined system with data 
determined by the substitute model. On the basis of the 
25 numerical value for the quality, the substitute model 
is adapted to be of as high a quality as possible. The 
substitute model adapted with regard to its quality is 
used for designing the technical system. 

This arrangement is particularly suitable for 
30 carrying out the method according to the invention or 
one of its developments explained above. 

Exemplary embodiments of the invention are 
presented and explained below with reference to the 
drawing, 

35 
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in which: 

figure 1 shows a block diagram which contains 
steps of a method for designing a technical system; 

figure 2 shows a schematic diagram of a 
5 recovery boiler; 

figures 3-5 show input variables, manipulated 
variables and output variables of the recovery boiler. 

Represented in figure 1 is a block diagram 
which contains steps of a method for designing a 

10 technical system. In a step 101, a substitute model is 
formed on the basis of measurement data. This 
substitute model is preferably a regression model. To 
adapt the substitute model created in step 101 to the 
measurement data, that is to perform a refinement of 

15 the substitute model, so that the measurement data 
describe the substitute model in adequate 
approximation, a numerical value for the quality of the 
substitute model is determined in a step 102. This 
numerical value is determined by comparing the 

20 measurement data of the predetermined system with data 
determined by the substitute model. Each item of 
measurement data preferably receives a numerical value 
for the quality, which numerical value characterizes 
the deviation of the item of measurement data from the 

25 associated value determined by the substitute model. 
The sum of all the numerical values for the quality for 
all the measurement data determines an overall quality 
for the substitute model. In a step 103, the quality 
is maximized by minimizing the numerical value for the 

30 quality (or a negative quality for the coincidence of 
the substitute model with the predetermined system) . 
Once an appropriately high quality for the substitute 
model has been determined, this substitute model is 
used for designing the technical system in a step 104. 

35 The designing may constitute both a new design (cf. 
step 105) or an adaptation of an already existing 
technical system (cf. step 106). 
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Figure 2 shows a schematic diagram of a 
recovery boiler. An exemplary embodiment of the method 
described above is illustrated below on the basis of 
the example of a "recovery boiler". 
5 In the paper and pulp industry, various 

chemicals and also heat and electrical energy are 
required for the digestion of pulp. The recovery 
boiler can be used to recover the chemicals used and 
also thermal energy from a concentrated process liquor 

10 (black liquor) . A measure of the recovery of the 
chemicals is of decisive significance for the cost- 
effectiveness of the overall plant. 

The black liquor is combusted in a char bed 
201. As this happens, an alkali fusion is formed, 

15 flowing away via a line 202. In further process steps, 
the chemicals used are recovered from the constituents 
of the alkali fusion. Released heat of combustion is 
used for generating steam. The combustion of the waste 
liquor and consequently the recovery of the chemicals 

20 begins with the atomization of the black liquor via 
atomizer nozzles 204 into a combustion chamber 203. 
Particles of the atomized black liquor are dried as 
they fall through the hot flue gas. The dried liquor 
particles fall onto the char bed 201, with a first 

25 combustion and a chemical reduction taking place. 
Volatile constituents and reaction products pass into 
an oxidation zone, in which oxidizing reactions proceed 
and in which the combustion is completed. 

Important targets for the controlling of the 

30 recovery boiler are the steam production for obtaining 
energy, the maintaining of emission values from 
environmental aspects and the efficiency of the 
chemical reduction. 
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The combustion operation, and consequently the 
targets, are controlled in particular by the air supply 
on three levels (Primary Air (PA) , Secondary Air (SA) , 
Tertiary Air (TA) ) . The overall process is subject to 
5 numerous influences, which have to be taken into 
account in the modeling: 

a) the measurement of the variables is often 
subject to strong fluctuations; 

b) unmeasured and unmeasurable influencing 
variables exist; 

c) any alteration of the parameters which can be 
set causes transient phenomena; 

d) the technical plant becomes soiled and is 
cleaned at predetermined intervals, having the 
result each time of a temporal drift in the 
system behavior. 

The measured variables of the overall process 
2 0 are divided into input variables (cf. figure 3) and 
output variables (cf . figure 5) . Measured values are 
stored every minute. Four of the input variables are 
at the same time also manipulated variables (also: 
parameters which can be set; cf . figure 4) . The 
25 manipulated variables are to be regarded essentially as 
free parameters of the overall process which can be set 
independently of one another. Some of the other input 
variables are more or less dependent on the manipulated 
variables. According to one target, in the case of the 
30 recovery boiler the "BL Front Pressure" and "BL Back 
Pressure" variables are always to be controlled such 
that they are the same. The four manipulated variables 
(cf. figure 4) are preferably to be stored as 
manipulated variables (with the desired, preset value) 
35 and as input variables (with the measured, actual 
value) . 

In the case of the recovery boiler, one problem 
is that of meeting targets which are determined in 
dependence on the parameters to be set and 

40 
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are defined by means of measured variables. A three- 
stage procedure is chosen here to solve the problem: 



1. The targets to be considered are modeled by 
5 stochastic methods, these models being updated by new 

measurements (data-driven, empirical modeling) . For 
this it is advisable to use not just a single model 
but instead global models for the identification of 
areas of interest in a parameter space determined by 
10 the targets and local models for the exact 

calculation of optimum operating points. The models 
used are assessed by measures of quality. 



, If the models considered are not sufficiently 
accurate (measure of quality) on account of the data 
situation, new operating points are evaluated on a 
specific basis to improve the model (experimental 
design) . Furthermore, by the use of global 

stochastic optimizing methods with respect to the 
targets, attractive areas are identified in 
dependence on the current global model. 

. For the local optimization, local models are devised 
and, if appropriate, the available sets of data are 
extended on a specific basis (experimental design) . 



The targets constitute physical-technical or 
business-management criteria, which generally have to 
conform to boundary conditions and/or safety 

30 conditions. Often a number of these criteria have to 
be considered simultaneously. The use of a stochastic 
model can be used in particular to simulate in the 
computer the target variables to be optimized and their 
dependence on the parameters to be set. This is 

35 necessary whenever measurements are very cost-intensive 
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or very time-consuming. In the case of safety 

requirements, possible hazardous situations can be 
avoided. 

In the case of the recovery boiler, an online 
5 optimization which is based on a number of data is 
necessary because the physical-chemical processes 
cannot be quantitatively modeled with sufficient 
accuracy and because the behavior of the plant is 
subject to fluctuations in the course of operation. 
10 Knowledge of this behavior must be continually extended 
by selective choice of new operating points. 
Therefore, the already described three-stage procedure 
of stochastic modeling and mathematical optimization is 
recommendable as part of online optimization. 

15 

DESCRIPTION OF THE INPUT VARIABLES 

The a input variables (a e N, N: set of 
natural numbers) are generally dependent on n actuated 
20 variables n e N and on random effects. They can be 
described as follows: 

Let (Q,S, P) be a probability space and be a 
Borel a algebra over /f^ (J2'f' set of real numbers) for 
each V e N. The input variables are represented by 
25 means of a B"^ x S - B^-measurable mapping cp: 

cp : R" X Q -> ^ 

The domain of the mapping cp is a Cartesian 
30 product of two sets. If the respective projections 
onto the individual sets are considered, the following 
mappings are obtained: 



cpx^ ^ -> R^, 0) <p(x, 0)) for all x e R"^ (2), 

35 
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<p^:R^ rS, X (p(x,co) for all co e Q 



(3) . 



|(px'" X e -'-^ ^ stochastic process with an index set 

and a mapping cp" is, for each event co e Q, 
a path of this stochastic process. 

In the case of the recovery boiler, n=4 and 
a=14 (after elimination of the "BL Back Pressure" 
variable) . 

On account of the required measurability of the 
mapping cp^, for each x e the mapping 9^ is a random 
variable. Under suitable additional preconditions, 
expectation values and higher moments can be 
considered. This approach permits the step from 

stochastic models to deterministic optimization 
problems. In the case of a deterministic optimization 
problem, the target function can be set directly by 
means of a variable, whereas the stochastic variable 
influences the target function but does not permit a 
specific setting. 



2 0 DESCRIPTION OF THE OUTPUT VARIABLES 



The process model M of the recovery boiler is 
described as a function in dependence on the input 
variables and further random effects. In this case let 
25 (Q,S, P) be the above probability space. The process 
model M is then a x S - B^-measurable mapping: 



30 where b denotes the number of output variables. 

Since the recovery boiler is subject to a 
cyclical temporal drift (from cleaning phase to 
cleaning phase) , a description with a time parameter is 
also conceivable. The 



35 
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output variables can be represented by B'' x S - B"^- 
measurable mappings 

Vj/ : R" X n ^ 
(x, ©) m(((>(x, w), (d) _ 

If the respective projections onto the individual 
sets of the domain are considered, the following 
mappings are obtained 



Vij,: Q R^, ® v(;(x, co) for all x e R" 
V":R" R^ X ^ ^;(x,(o) for all cd e Q 



; X € R'^I "'"^ ^ stochastic process with an index set 

%: ^ 

and the mapping xj/™ is, for each oo e Q, a path of 
this stochastic process. 

In the case of the recovery boiler, b=15. 
15 The fact that in the definition of v|/ a 

distinction is not drawn between the events oo used does 
not mean that there is any restriction, since Q can be 
represented as a Cartesian product of an Di and an . 
The above representation consequently also comprises 
20 the model: 

\{/ : R" X X Q2 Rt' 

(x, coi, C02) m((p{x, a>i), 032) 



(9) , 

(10) . 



DESCRIPTION OF THE AVAILABLE SETS OF DATA 

25 
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With the descriptions in the two foregoing 
sections, the input variables and the output variables 
can be combined together to form measured variables (= 
measurement data) O. O is a B" x S - B^'-measurable 
5 mapping with m = a + b and 



I ^ ( <p(x, (£>)] 

If the respective projections onto the 
individual sets of the domain are again considered, the 
10 following mappings are obtained: 



Ox: Q R^, CO -» (t>(x, 0)) for all x e (13), 

<I>*":R'' a)(x,©) forall CO € Q (14). 



|<I>x; X € R"| "'"^ ^ stochastic process with an index set 

X", and the mapping O'" is, for each cd e Q, a path of 
15 this stochastic process. 

For each chosen related variable tupel x, in 
the case of the recovery boiler many realizations of 
are determined and stored, i.e. for each Xj e 
numerous realizations 

20 

Ojk := o(xj,©jk) (15) 
with CD j k e Q ; k = 1,2, . . . , V j ; 



V j e N ; j = 1,2, . . . , u ; u e N 
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are considered. The stored sets of data Djk of the 
recovery boiler are consequently (n+m) tupels : 

^jk = [<^-y^ ' ^ = lr2 Vj ; j = 1,2 u (16) . 

5 

In this case, Djiki is stored before Dj2k2 if 

Gi < 32) V (Oi = 32) A (ki < k2)) 

applies . 

10 DATA COMPRESSION BY CLASSIFICATION OF THE PARAMETERS 

Since for each manipulated variable tupel x there are 
generally a number of realizations of O^, a 
classification of the parameters by forming arithmetic 

15 mean values is a suitable operation as a first step of 
the statistical data analysis on account of the complex 
stochastic properties of the process to be considered. 
Moreover, sets of data which are obviously erroneous 
are picked out. A set of data which is obviously 

20 erroneous is, for example, a physically impossible 
measurement which cannot occur in reality, in 
particular on account of a setting made. 



Procedure : 

25 

1. Sets of data in which the "BL Front Pressure" 

variable is not equal to the "BL Back Pressure" 

variable are picked out, since these two values must 
be equal according to the default of the plant 

30 control. The data loss is very small. 



2. The sets of data are divided into classes in which 

the four setting parameters (PA, SA, TA, BL Front 

Pressure, see above) are constant in temporal 
35 succession. 
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i.e. the jth class comprises the data sets Dj^. 

3. Classes in which there are fewer than 30 data sets 
are picked out in order that transient phenomena have 
no great influence. 

4. For each class, an arithmetic mean value Oj and an 
empirical standard deviation Sj are determined for 
all the measured variables: 



1 



5. Classes in which the mean values for the variables 
PA, SA, TA or BL Front Pressure are too far away from 
15 the corresponding setting parameters are picked out. 

Therefore, in these classes it was not possible to 
reach the setting values. 
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STATISTICAL CHARACTERISTIC VARIABLES FOR THE GIVEN 
CLASSES AND THEIR GRAPHIC REPRESENTATION 

Apart from the arithmetic mean values and the 
5 empirical standard deviations which were determined for 
the individual classes, a common standard deviation s 
is also determined according to 



205) and v stands for the sum of Vj, i.e. v is the 
number of all the measured values used (here 38 915) . 

15 LINEAR REGRESSION MODELS FOR FUNCTION APPROXIMATIONS 

For each measured variable (item of measurement 
data) O^"-' (i=l, 2, . . . ,m) , a linear regression model is 
calculated on the basis of the arithmetic mean over the 
20 classes in dependence on the quadratic combination of 
the four setting parameters. In the following 

representation, x e , where 




(19) 



10 



where u stands for the n^umber of classes (here 



25 




Primary Air (PA) 
Secondary Air (SA) 
Tertiary Air (TA) 

Black Liquor (BL) Front Pressure 



(3) 



5 
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applies. u e N denotes the number of classes. 
Each measured variable O*^' is modeled by 



O^^^x, to) = ajr(x) + ei(to) 



with ai € R^^. Here the following applies: 



(21) 



fel^ ^2r ^3r ( 1, ^1, ^2. ^3' ^1 ^2' Ci' 

;iC2.- ^1^3. ^lC4. ?2;3. <;2;4' ^3^4 ^ ^ (22 ) , 

i.e. polynomials of the second degree are adapted to 
10 the measurement data, and 

ei : Q R (23) 



is a random variable with the expected value 0. 

The vector ai is determined by the method of 
least squares, but the arithmetic means { 



are used instead of the original data se 



This procedure is suitable, since linear regression 
models estimate in particular expected values. This 
results in the following minimization problem: 



^1 

.(15) 



20 
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Let ai be the optimum point of the quadratic 
minimization problem from equation (24) . Furthermore, 
the following applies: 



To validate the regression theorem, a 
coefficient of determination is calculated according 



{_ ^ YiYi - uyf 



B(i) 



The closer R i is to 1, the better the dependent 
variable is represented by the independent variables 
15 (o < Ri < l) 

In addition, a maximum E'^'max for an absolute 
value of the deviation of the data from the model is 
specified as 



E&L:= max |#^-yW|| 
j = l,..-,u 11 IJ 
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E 



(i) 



90% is that value below which at least 90% of 



the absolute values of the deviations of the data from 
the model lie. By analogy with this, E*^'eo% is that 
value below which at least 8 0% of the absolute values 
of the deviations of the data from the model lie. With 
the optimum point ~al of the minimization problem 
according to equation (24), a model O*"-' of the expected 
value of the measured variable O*^' can be specified as 



In particular, the gradient VO*^' can be 
analytically specified by 



^(i) := ^ R 



(30) , 




(31) . 



va^^^(x) = ^ (x) ■ ai for all 



X € R' 



,n 



(32) . 
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1- A method for designing a technical system, 

5 a) in which measurement data of a predetermined 

system are described on the basis of a 
substitute model; 



b) in which a numerical value for the quality of 
the substitute model is determined by comparing 
the measurement data of the predetermined 
system with data determined by the substitute 
model; 



c) in which the substitute model is adapted from 
the numerical value for the quality to be of as 
high a quality as possible; 

d) in which the substitute model adapted with 
regard to its quality is used for designing the 
technical system. 

2. The method as claimed in claim 1, in which the 
substitute model is a regression model. 

3. The method as claimed in claim 1 or 2, in which 
the quality is determined on the basis of a mean square 
deviation of the measurement data from the data 
determined by the substitute model. 

4 . The method as claimed in one of the preceding 
claims, in which the measurement data are sorted 
according to their quality, with respect to the 
deviation of the latter from the data determined by the 
substitute model, and a predetermined number of n% of 
the worst measurement data are picked out. 

5. The method as claimed in claim 4, in which the 
n% of the worst measurement data are not picked out 
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if they lie in a continuous range. 

6 . The method as claimed in one of the preceding 
claims, in which the amount of measurement data is 
reduced in the course of a preprocessing operation. 

7. The method as claimed in claim 6, in which the 
preprocessing operation comprises a classification of 
the measurement data. 

8 . The method as claimed in one of the preceding 
claims, in which the data obtained by means of 
designing are used for controlling a technical plant. 

9. The method as claimed in claim 8, for the 
online adaptation of the control for the technical 
plant . 

10. An arrangement for designing a technical 
system, with a processor unit which is set up in such a 
way that 

a) measurement data of a predetermined system are 
described on the basis of a substitute model; 

b) a numerical value for the quality of the 
substitute model is determined by comparing the 
measurement data of the predetermined system 
with data determined by the substitute model; 

c) the substitute model is adapted from the 
numerical value for the quality to be of as 
high a quality as possible; 

d) the substitute model adapted with regard to its 
quality is used for designing the technical 
system. 
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Abstract 

Method and arrangement for designing a technical system 

To achieve the object, a method for designing a 
technical system in which measurement data of a 
predetermined system are described on the basis of a 
substitute model is specified. A numerical value for 
the quality of the substitute model is determined by 
comparing the measurement data of the predetermined 
system with data determined by the substitute model. 
On the basis of the numerical value for the quality, 
the substitute model is adapted to be of as high a 
quality as possible. The substitute model adapted with 
regard to its quality is used for designing the 
technical system. 
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